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We consider the stability of a compressible shear flow separating two streams of different speeds and temperatures. The 
velocity and temperature profiles in this mixing layer are hyperbolic tangents. The normal mode analysis of the flow 
stability reduces to an eigenvalue problem for the pressure perturbation. We briefly describe the numerical method that 
we used to solve this problem. Then, we introduce the notions of the absolute and convective instabilities and examine the 
effects of Mach number, and the velocity and temperature ratios of each stream on the transition between convective and 
absolute instabilities. Finally, we discuss the implication of the results presented in this paper for the heliopause stability. 
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1 Introduction 

The theory of stabiHty of shear flows has been the subject 
of an intense study for the past century due to its numerous 
applications in engineering, geophysics, and astrophysics. 
Velocity shear generates a dynamical instability, and the 
most common and fastest one is the Kelvin-Helmholtz (KH) 
instability which occurs in fluid and plasma inviscid shear 
layers. It intervenes in many phenomena in fluid dynamics 
(Blumen et al. 1975; Jackson & Grosch 1989; CaiUol 2005), 
in space, astrophysical and laboratory plasmas (Miura 1982; 
Jones et al. 1997; Mifls et al. 2000; Ruderman 2000; Taro- 
yan & Erdelyi 2002, Hasegawa 2004) when steep velocity 
gradients emerge. The most relevant examples are the stabil- 
ity of the interface between the solar wind and the planetary 
magnetospheres, the interaction between adjacent streams 
of different velocities in the solar wind and polar cusps, the 
dynamic structure of cometary tails, and the interaction be- 
tween the solar wind and the interstellar medium. 

In space and astrophysical plasmas, the total velocity 
jump across the velocity shear layer can have large sonic 
Mach number, so that taking into account the compressibil- 
ity of plasmas is essential. KH instability weakens as the 
convective Mach number increases (Mach number in the 
frame moving with a perturbation), while oblique distur- 
bances become more and more unstable and dominant. 

In astrophysical applications, the account of magnetic 
field is very often important. In that case, the fluid motion is 
described by the magnetohydrodynamic (MHD) rather than 
hydrodynamic (HD) equations. Studying the MHD KH in- 
stability is much more complicated than the HD one due to 
the coupling between KH instabilities and Alfven field line 
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resonances (Hollweg et al. 1990; Yang & Hollweg 1991; 
Ruderman & Wright 1998; Taroyan & Erdelyi 2003). In fliis 
paper, we study the HD KH instability postponing the mag- 
netic field account for future work. 

The study of stability of compressible hydrodynamic 
shear flows has a long history. Landau (1944) analytically 
studied the discontinuous flow generated by a vortex sheet 
and proved its stability when the Mach number M is greater 
than \/2 with respect to two-dimensional disturbances. Blu- 
men et al. (1975) showed that the inviscid flow with the ve- 
locity profile U{y) = tanh(y), where y is the Cartesian co- 
ordinate perpendicular to the flow direction, is temporally 
unstable for any M. 

The normal mode analysis used in these studies revealed 
many important properties of the KH instability of shear 
flows. However, it failed to predict if a finite portion of flow 
looks stable or unstable in a fixed reference frame. The rea- 
son is that the normal mode analysis deals with periodic 
perturbations infinite in at least one spatial direction, while 
real perturbations always occupy a finite spatial domain. To 
answer this question, we need to solve the initial value prob- 
lem for an arbitrary perturbation bounded in a finite domain. 
Then two scenarios are possible. In the first scenario, the 
perturbation exponentially grows with time at any fixed spa- 
tial position. This type of instability is called absolute. In the 
second scenario, the perturbation exponentially grows with 
time, but in the same time it is convected out of the observa- 
tional domain so fast that eventually it decays at any fixed 
spatial position. This type of instability is called convective. 

The notions of the absolute and convective instabilities 
were first introduced in plasma physics (see Briggs 1964; 
Bers 1973). Later, it started to be used in hydrodynamics 
(e.g., Kulikovskii & Shikina 1977; HueiTe & Monkewitz 
1985; Brevdo 1988). Initially, the absolute and convective 
instabilities were studied for shear flows of incompressible 
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fluids. Pavithran & Redekopp (1989) analysed the convec- 
tive and absolute instabilities in a shear flow of a compress- 
ible fluid when the velocity and temperature have hyperbo- 
lic-tangent profiles. Jackson & Grosch (1990) carried out a 
similar study for a boundary layer flow where the velocity 
field is described by the hyperbolic tangent, and the tem- 
perature is linked to the velocity by the so-called Crocco 
relation. In both papers, the analysis is restricted to the sub- 
sonic regime. Terra-Homen & Erdelyi (2003) determined 
the transition from the convective to the absolute instability 
in a shear flow of a viscous compressible fluid by solving 
the initial value problem for the full nonlinear system of 
hydrodynamic equations. In their study, they considered a 
shear flow with a constant temperature and the hyperbolic- 
tangent velocity profile. They also restricted their analysis 
to subsonic flows, and took the Reynold number equal to 
lO'^. With these constraints, the stabilizing effect of viscos- 
ity was stronger than that of compressibility. 

The aim of this paper is to study the transition from the 
convective to the absolute instability in the same shear flow 
as one considered by Pavithran & Redekopp (1989), how- 
ever with the Mach number attaining supersonic values. The 
paper is organized as follows. In the next section, we for- 
mulate the problem, and derive the governing equation and 
the boundary conditions for the pressure perturbation. In 
Sect. [3] we describe the spectral method used to solve the 
eigenvalue problem determining the flow stability with re- 
spect to normal modes. In Sect. HI we present the results of 
the study of the flow absolute and convective instabilities. 
In Sect.|5] we discuss the implication of the results obtained 
in this paper for the heliopause stability. Sect.|6]displays the 
summary of the results and our conclusions. 

2 Formulation 

We consider the instability of a two-dimensional and un- 
bounded compressible ideal shear flow. The current is in the 
cc-direction. Both velocity and temperature magnitudes are 
functions of y in Cartesian coordinates x,y,z and are given 
by 

U = i[l + /3c/ + (l-/3t/)tanli(r;)], (1) 

T = ^[l + f^T + {I ~ (3t) tanh(ry)] , (2) 

where ?/ is the cross-stream variable in the Howarth-Dorod- 
nitzyn transformation, 77 = J^^ 1/T{y') dy' (Jackson & 
Grosch 1989). The velocity and temperature are made di- 
mensionless using their values in the fast stream (y 00). 
The spatial variables and time are also dimensionless and 
measured in units of the half-width of the shear layer yo, 
and yo divided by the fast stream velocity, respectively. The 
quantity (3u is the ratio of the speed of the slow stream to 
that of the fast stream, and /3t is the same but for tempera- 
ture. We always have 



-l</3c/<l, /3t>0, 
U^Pu, T^Pt at ?7 ^ - 




Fig. 1 Profile U{y) for various values of f3u, /3t = 2 
(solid line) and (3t — 10 (dotted fine). 



Fig. 2 Profile T{y) for /3t = 0.1, 0.5, 1, 2, 5 and 10. 



The velocity and temperature profiles are shown for differ- 
ent values of Pu and Pt in Fig. [T] and Fig. |2] respectively. 
Flows of this type represent transition layers between two 
streams with constant parameters. The KH instability re- 
sults in mixing of the fluids from both streams inside the 
transition layer. This is why such flows are called mixing 
layers. 

The flow field is perturbed by introducing wave distur- 
bances whose amplitudes are functions of for example, 
the pressure perturbation is 

p = 11(77) eyip[i{ax + Pz ~ (lit)] , 

with a and P the wavenumbers and d) the frequency. The 
equation governing 11 (Jackson & Grosh 1989) is 

2t/'n' 



u 



- T[{c 



P^)T 



a^M'^{U - cf]n = 0, 



(3) 



where c — oj/a, and the prime indicates the derivative with 
respect to rj. Equation (|3]l is easily turned into an analogous 
two-dimensional disturbance equation by Squire transform 



a 



-P' 



a = a cos t 



P = a sin ( 



aM = aM , c = c , ili — ui cos 9 , 

where 9 is the angle between the disturbance wavenumber 
and the mean flow direction. Applying this transform to ^ 
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yields 



U 



- a^T[T - (U - cY]n = . 



(4) 



Whatever n(Q;, c) is a solution of (|4|i, n(— a, c) and 
Il*{a*,c*) are also solutions, so that we can restrict the 
analysis in the (a, c)-plain to a,. > and q > 0, where 
the indices r and i indicate the real and imaginary parts of a 
quantity. 

The boundary conditions for 11 are obtained by consid- 
ering the limiting forms of ^ as rj ^ ±oo, that is by taking 
[/' = in (IDi. The required solutions are exponentials of the 
form 

n exp(=Fri±?;) at rj ^ ±oo , (5) 

where n\ = a^[l - M^{1 - cf] , 

n't = a^p^lpT-M^iPu-cf]. 

The linear theory leads to consider small-amplitude pertur- 
bations, which are weakly unstable oscillations character- 
izing sonic waves escaping from the shear layer, in which 
^{vgy) > when y > and ^{vgy) < when y < 0; 
Vgy is the y-component of the group velocity. This warrants 
an exponential decrease of the wave energy away from the 
shear layer The non-radiation condition for unstable oscil- 
lations comes down to "31(^1^) > and 3fi(il_) > 0, where 
3? indicates the real part of a quantity. A neutral mode with 
a subsonic convective Mach number decays exponentially 
with 1] but radiates outward with a supersonic convective 
Mach number like a sound wave. An unstable supersonic 
mode may exponentially decay but more slowly than the 
subsonic modes. 

As the profiles {U, T) have constant "tails," the unbo- 
unded-flow problem can come down to a finite-boundary 
problem in the 77-range [?7_ , 77+] such as — J7(±oo)| <C 1 
as I77I > min[— 77_, 77+] (Keller's boundaries). Conditions 
Q therefore become 



3 Spectral Method 



dr] 



0. (6) 



When c is given, equation (HJi together with the boundary 
conditions (|6]l constitute an eigenvalue problem, a being the 
spectral parameter. The solution of this problem gives a as 
a function of c, or, assuming that this function can be in- 
verted, c and, consequently, w as a function of a. To solve 
this problem we use the following numerical method. First 
we introduce a new cross-stream variable Z such as 



-l< Z <l 



(7) 



where -q^ = —rj+ and < /i < 1. Equation ^ is rewritten 
as 



fa 

dZ2 



3(1-^);^ + 



{ai^+EfT[T 



dU/dZ 

U-c 



dn 

dz 

\2 



where E ^ fi + 3{1 - fj.)Z'^ , 

^^i(l-/3y)r,+i?[l-tanh(,7)2] 

The boundary conditions (|6]l take the form 

— ±a'q,Ew±U = at Z = ±1, 
dZ 



(9) 



(10) 



where zu± — il±/a. The problem comes down to the 
search for the eigenvalues A = ar]+E{l) > 0) 

and eigenfunctions 11, the wave speed c being given. An 
approximate solution is obtained by expanding 11 in a fi- 
nite Chebyshev series. The Chebyshev collocation method 
is used to discretize the problem. This method was tested 
by a direct integration of (|4|i using the fourth-order Runge- 
Kutta scheme in the same way as in Blumen (1970). 

4 Convective and Absolute Instabilities 

The general evolution of an initial perturbation with time 
can be formally expressed by a Laplace-Fourier integral 



q{x,y,t) 



ia-\-OQ 

e-'^' du 

a — 00 



-00 D{a,uj) 



e'^'^da, (11) 



where a is defined to be larger than the largest growth rate 
of any mode for real a, and q is the perturbation of any 
quantity. The dispersion relation D{a, uj) and the function 
T{a,uj,y) are analytic functions of a and w. The function 
T{a, uj, y) depends on the initial perturbation. It is not im- 
portant for what follows because the asymptotic behaviour 
of q{x,y,t) as i 00 is completely determined by 
D{a,uj). Note that the initial-value problem is well-posed 
because the growth rate of unstable modes in an inviscid 
compressible shear flow is bounded by (Chimonas 1970) 



(12) 



To study the asymptotics of q{x, y, t) , we used the Briggs' 
method (Briggs 1964). In accordance with this method, this 
asymptotics is closely related to the double roots of D{a, uj) 
considered as a function of a. The double a-roots satisfy the 
system of equations 



D{a,u;) = 0, 



dD 
da 



(a, uj) = Q . 



(13) 



c)1n = 0, 



(8) 



Due to the restriction on the paper length, we do not de- 
scribe the details of solutions (a, lo) of this system of equa- 
tions. They will be given elsewhere. In this paper, we will 
only outline the main results that we obtained. 

Our main results are presented now. The figure[3]shows 
the boundaries between the regions of parameters corre- 
sponding to the absolute and convective instabilities in the 
(M, /?(7)-plane for different values of /3t. A general trend 
is apparent whatever (3t- There exists a critical Mach num- 
ber Afcr whose value depends on (3t such as when M < 
M„, increasing the Mach number increases the magnitude 
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Fig. 3 The boundaries between the regions of parameters 
corresponding to the absolute and convective instabilities in 
the (M, /Jjy) -plane for different values of (3t. 'AU' stands 
for the absolute instability, and 'CU' for convective insta- 
bihty. 

of counterflow necessary to cause the absolute instabil- 
ity. But as AI > M„, the situation is opposite: the larger 
the Mach number is, the smaller the magnitude of counter- 
flow is. Cooling the slow stream, i.e. decreasing Pt, extends 
the domain of the absolute instability. For (3t sufficiently 
small, the shear layer may become absolutely unstable in 
the subsonic range even when the streams are co-flowing, 
i.e. f3u > 0. In particular, the threshold value of Pu be- 
comes positive for /3t — 0.1 when AI < 0.8, and reaches 
a value of order 10^"^. A finer accuracy and a stronger reso- 
lution of the spectral method were used for these very small 
values of (3u. All curves tend to zero as Af ^ oo. This 
implies that the flow is absolute unstable for any value of 
(3u < if AI is sufficiently large. Figure|4]shows the transi- 
tion frequency as a function of AI for chosen values of 
A similar separation of the evolution of ujr with AI occurs 
according to whether AI is less or greater than Alcr. In the 
first range, uJr is quasi-constant, then it decays rapidly for 
M ^ Mcr. When the motion is supersonic, the decrease 
is smoother. When Pt — 0.1 however, the behaviour of 
UJr is different, it is increasing up to M = Mcr and then 
decreasing. The hotter the slow stream is, the smaller the 
transition frequency is. The asymptotic frequency for infi- 
nite Mach number is zero. Figure |5] displays the evolution 
of the wavenumber as a function of AI. This does not vary 
substantially for AI < Mcr when (3t is of order one, but 
it increases steeply when Pt is small. The wavenumber 
increases as /3t decreases. We expect the high instability of 
the cooled-slow-stream mixing layers to be attenuated by 
viscous and thermal diffusions. 



5 Application to Heliopause Stability 

Collision of the solar wind with the interstellar medium re- 
sults in the interaction region consisting of two shocks and a 
tangential discontinuity between them (Baranov et al. 1970, 
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Fig. 4 Absolute/convective instabilities transition in the 
(w,-, AI) plane for various values of Pt- 
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Fig. 5 Absolute/convective instabilities transition in the 
(ttr, AI) plane for various values of Pt. 

1976, 1979; Baranov 1990). This region is shown in Fig.|6l 
The external shock, called the bow shock, decelerates the 
supersonic flow of the interstellar medium. The internal 
shock, called the termination shock, decelerates the super- 
sonic flow of the solar wind. The tangential discontinuity, 
called the heliopause, separates the decelerated solar wind 
and intersteUar medium flows. 

The heliopause stability was a subject of an intensive 
study for more than two decades. In the early studies, the 
normal mode analysis combined with the local analysis (the 
WKB approximation) was used (e.g., Fahr et al. 1986; Bara- 
nov et al. 1992; Chalov, S.V. 1996). Recendy, Ruderman et 
al. (2004) made an attempt to study the absolute and con- 
vective instability of the heliopause. However, since these 
authors used the approximation of an incompressible fluid, 
they only managed to study the stability of near flanks of 
the heliopause (the region corresponding to 16*1 < 30° in 
Fig.|6]l. At the near flanks, the plasma flows at both sides of 
the heliopause are strongly subsonic, so that the approxima- 
tion of an incompressible fluid can be used to describe the 
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Fig. 6 Schematic of the heliosphere 

plasma motion. However the flow Mach number monoton- 
ically increases with the distance from the apex point, and 
eventuaUy the flow becomes supersonic at far flanks, which 
implies that the account of compressibility is of crucial im- 
portance. An immediate implication of the results presented 
in the previous section is that the whole heliopause is only 
convectively unstable, at least with respect to short-wave- 
length perturbations when the local analysis is valid. Really, 
even when the ratio of temperatures at the two sides of the 
heUopause, /3t, is very small (which is usually the case), 
the instability is absolute only if Pu 10^''. However, in 
accordance with all numerical models of the solar wind- 
intersteUar medium interaction 10~^ (e.g. Baranov 

et al. 1979). 

6 Summary and Conclusions 

In this paper, we have studied the absolute and convective 
instabilities of a compressible shear flow with the hyper- 
bolic-tangent velocity and temperature profiles. The stabil- 
ity properties of the flow are determined by three dimen- 
sionless parameters: the Mach number M = M cos 9 calcu- 
lated using the velocity and sound speed of the fast stream, 
and the ratios of the velocities f3u, and temperatures /3t, 
of the fast and slow streams. We plotted the boundaries be- 
tween the regions of parameters corresponding to the abso- 
lute and convective instabilities for different values of /3t- 

Our main results are the following. Each value of /3t 
imposes a threshold value of the Mach number Ma. When 
M < Mcr, increasing the Mach number increases the mag- 
nitude of counterflow (i.e. increases the value of — /Jy) nec- 
essary to cause the absolute instabihty. This result was pre- 
viously obtained by Pavithran & Redekopp (1989) and Jack- 
son & Grosch (1990) for subsonic flows. When M > Mcr, 
increasing the Mach number decreases the magnitude of co- 
unterflow necessary to cause the absolute instability. Co- 
oling the slow stream, i.e. decreasing (3t, extends the do- 
main of the absolute instability. For very small values of Pt 
the absolute instabihty is even possible when the streams 
are co-flowing, i.e. pu > 0, however only for very small 



values of /3u, j3u < 10~^. The magnitude of the counter- 
flow necessary to cause the absolute instability tends to zero 
i-Pu 0) when M ^ oo. 

We discussed the implication of the results obtained in 
this paper on the stability of the heliopause, which is the 
surface separating the flows of the intersteUar medium and 
the solar wind plasma in the solar wind-interstellar medium 
interaction region. Our conclusion is that the heliopause is 
only convectively unstable, at least with respect to short- 
wavelength perturbations. 
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